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We introduce a conceptually simple and experimentally feasible method to realize and detect 
photonic topological Chern insulators with one-dimensional circiut quantum electrodynamics lattice. 

By periodically modulating the couplings in this lattice, we show that this one-dimensional model 
can be mapped into a two-dimensional Chern insulator model. In addition to allowing the study of 
photonic Chern insulators, this approach also provides a natural platform to realize experimentally 
Laughlin’s pumping argument. Remarkably, based on scattering theory of topological insulators 
and input-output formalism, we hud that both photonic edge state and topological invariant can be 
unambiguously probed with a simple dissipative few-resonator circuit-QED network. 

PACS numbers: 85.25.Am, 85.25.Cp, 03.65.Vf 


Introduction. The recent rapid experimental develop¬ 
ments in circuit quantum electrodynamics (circuit QED) 
have turned this system into one of leading platforms 
for studying quantum optics and quantum computation 
[I[. This system possesses high coherence superconduct¬ 
ing qubits [ll and well controllable coupling parameters 
@ 1 13- Putting qubits and microwave resonators into a 
lattice further allows this system to be used for solid state 
quantum simulation. Quantum state can be easily ma¬ 
nipulated and detected at single-site level in such lattice. 
Combined with on-site nonlinearity or photon blockade, 
circiut QED lattices have been widely explored as quan¬ 
tum simulators in the past years for invest^ating photon- 
or polariton-based many-body physics [a-lTj . Experimen¬ 
tally, mimicking quantum spin models with (mbit arrays 
recently has been successfully demonstrated Q. 

Topological photonics nowadays has become a very ac¬ 
tive area of research Q. Photonic topological insulator 
was firstly predicted in a two-dimensional photonic crys¬ 
tal and subsequently has been extensively studied 
[lll - ll5l| . Based on engineering artifical magnetic field and 
spin-orbit coupling, photonic quantum integer and spin 
Hall states have also been studied in two dimensional cou¬ 
pled resonators (l6l - l^ and linear circuit lattice [l^, [13 . 
Experimentally, the propagation of photonic chiral edge 
states has been observed [ll, [H, [13 ill, Ull . It would 
be highly desirable, in these experiments, to measure 
the photonic topological invariants. Due to the differ- 
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ent statistics of the carrier, photonic topological invari¬ 
ants cannot be measured using the same methods devised 
for electronic systems. Recently, how to detect photonic 
topological invariants has been studied based on prob¬ 
ing the Berry curvature [13l and the dynamics of edge 
states [2^,[23. However, although both circuit QED and 
topological photonics have been rapidly developed, the 
connection of them is less studied. 
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FIG. 1: (a) Setup for the one-dimensional circuit QED lattice, 
where the cavity input-output process has been explored to 
probe the edge state and topological invariant, (b) An artih- 
cial torus is constructed by combining the parameter 6 with 
the quasimomentum of the photonic lattice. 

In this paper, we propose a conceptually simple 
method to simulate a two-dimensional photonic topolog¬ 
ical Chern insulator with a one-dimensional transmission 
line resonator photonic lattice. By identifying a peri¬ 
odic parameter introduced in the system as the quasi¬ 
momentum in a second artificial dimension, we show 
that this one-dimensional photonic lattice can be mapped 
onto a two-dimensional lattice exhibiting the Chern in¬ 
sulator state. Compared to previous methods of engi- 
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neering photonic topological insulators, our method em¬ 
ploys a one-dimensional photon lattice, and it does not 
require mimicking gauge fields and spin-orbit coupling, 
which opens a new and simple route to study and probe 
high dimensional photonic topological states. It is also 
interesting to note that such one-dimensional system can 
provide a natural platform to realize Laughlin’s pumping 
argument [33 and to test scattering theory of topological 
insulators [^. Based on such feature, we generalize such 
scattering theory to take into account dissipation and 
demonstrate that both photonic edge states and topo¬ 
logical invariant still can be clearly probed from the final 
steady state. We believe that our results could simu¬ 
late more studies on the dissipation effect in topological 
states. It is also noteworthy to stress that the topolog¬ 
ical features emerge even in a dissipative few-resonator 
circuit-QED network. Topological invariant defined on 
a two-dimensional topological lattice can be unambigu¬ 
ously detected with a small-size one-dimensional simula¬ 
tor. 

One-dimensional circuit QED lattice. We start with 
a transmission line resonator lattice, as shown in Fig. 
1(a). The photon hopping between nearest neighbour 
resonators is mediated through the coupling capacitors 
and the connected flux qubits. Note that this cou¬ 
pling method has recently been demonstrated in the two- 
resonator circuit QED experiments Si- In both cases, 
the coupling has been designed to make sure that the res¬ 
onator lattice has the alternating hopping configuration. 
Each unit cell of this lattice has two resonators, labelled 
by a and b. The capacitively coupled resonator lattice is 
described by the Hamiltonian 


Hq ='^Jiai^bnJ2a:j^bn-ih.C; ( 1 ) 


where Ji and J 2 are the intra- and inter-cell hopping 
rates. 

For the qubit-assisted hopping, we assume that the two 
resonators within the same unit cell are both coupled to 
the flux qubit Qi, while the two resonators belonging 
to the two nearest-neighbour unit cells are both coupled 
with the flux qubit Q 2 - The purpose of this coupling 
is to provide an alternating parametric modulation on 
the hopping rates and the on-site energies. In the dis¬ 
persive regime, when all the qubits are in the ground 
state, the coupling between the resonator and the qubit 
can be removed, leading to an effective transmission res¬ 
onator lattice with photon hopping assisted by the con¬ 
nected qubits. Combined with the previous capacitively 
coupled resonator lattice, the total Hamiltonian of this 
cicuit-QED lattice (in a rotating frame with respect to 
the external driving frequency uJd and also in the interac¬ 
tion picture with respect to the qubit energy a;i^ 2 ) takes 


the form 

H = ^(Tl — + {J2 + + h.C 

n 

2 2 

+ A~ + Ac(a)(^a„ -f 

( 2 ) 

where gi and 32 describe the coupling strengths between 
the qubit Qi{Q 2 ) and the resonators a„ and bn, {bn and 
a„+i ), A = uii — uJd = cjd — UJ 2 is the detuning of the 
qubit energies, and Ac = uJc — uJd is the detuning of the 
resonator frequency. The qubit-assisted hopping and on¬ 
site modulation terms are introduced in order to map into 
the second dimension for the construction of the photonic 
Chern insulator. 

Two dimensional lattice mapping. To simulate the two- 
dimensional Chern insulator Hamiltonian [s^, we write 
the qubit-resonator coupling strengths in the above lat¬ 
tice Hamiltonian in a parameter space as 

91 = 90 sin(6»/2), 32 = 9o cos(6</2), (3) 

where the mixing angle 9 = 2 arctan( 5 i/g 2 ) and go = 
\Jg\ + gf. The parameter 9 is determined by the ratio 
between the coupling strength gi and 32 - Note that the 
coupling strengths between the flux qubit and the res¬ 
onators can be individually controlled through using su¬ 
perconducting quantum interferences (SQUIDs) devices 
and changing the external magnetic fluxes applied on the 
SQUIDs loops [ 13 . Then 9 can be engineered from 0 
to 2tt for subsequent two dimensional mapping. More¬ 
over, the topological feature demonstrated below in this 
model endows this system with topological protection, 
which allows our methods to be robust to practical de¬ 
formations in the parameters engineering. By substitut¬ 
ing the above equation into the total lattice Hamiltonian 
and further writing it in momentum space, one can get 
H = h{k)Ck, where Ck = {ak, bkY. The momen¬ 

tum density has the following form 

h{k) = ho + hxCJx + hyffy -I- h^a^, (4) 

where ho = Ac and h = {hx,hy,hz} = 
{2 J cos{kx), 2(5 sin{kx) — Je sin(0) siii{kx),Je cos(0)} with 
J ={Ji+ T2)/2, 5 = {Ji- J2)/2 and Je = gHA. ax,y,z 
are the Pauli matrices spanned by and bk- Interest¬ 
ingly, by associating the mixing angle 9 with the quasi¬ 
momentum ky in the second spatial direction, one can 
find that the above one dimensional circuit QED lattice 
can be exactly mapped into a two dimensional Chern 
insulator Hamiltonian. As plotted in Fig. 1(b), the x 
direction quasimomentum kx and the mixing angle 9 can 
form a two dimensional Brillouin zone kx € f0,7r] and 
0 e [0, 27r], which can be rolled into a torus |3J-[^ 
for analyzing the underlying topology in the artificial two 
dimensional lattice. 

Photonic Chern insulator. The topological properties 
of the model introduced above are captured by the Chern 
number of the Bloch band and the edge state spectrum. 
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FIG. 2: Edge state spectrum for Photonic Chern insultaor 
with (a) Chern number C = 1 when J = 0 and (b) Chern 
number C = 0 when S — 0.6Je- For Chern insulator, there 
are two edge states at the in-gap energy denoted by the red 
dashed line. Inset: the density distribution of the two edge 
states. The other parameter are chosen as J = Je and the 
lattice size L = 10. 


By mapping the two-dimensional torus to a spherical sur¬ 
face, the Chern number of the occupied ground band 
can be expressed as C = ^ f f dkxdd(dk^h x clgh) • h, 
where the unit vector h = (hx, hy, hz)/\h\ with \h\ = 
+ hi- Through substituting hx,y,z into above 
formula and one can get the Chern number of the ground 
band as 


C = 


if - Je <25 < Je 
otherwise 


( 5 ) 


One can change the hopping difference 5 to engineer 
the photonic topological phase transition. It is also 
worth pointing out that, when the coupling strength 
(72 = ~JoCOs(0/2), the ground state can be prepared 
as a Chern insulator with C = — 1. According to the 
bulk-edge correspondence, the appearance of edge state 
is a hallmark of topological insulator. In Fig. 2(a), we 
have plotted the edge state spectrum for the topological 
insulator. There is one pair of edge state at the in-gap 
energy. The density distribution for the left and right 
edge states have been plotted in the inset figure of Fig. 
2(a). One finds that there is one edge state localized at 
each edge. For the topological trivial insulator shown in 
Fig. 2(b), there is no edge state in the gap. 

Scattering formulation of topological invariant. Based 
on Laughlin’s pumping argument, recent scattering 
theory of topological insulators shows that topological 
invariant can be described by the reflection matrices at 
the Fermi level M- The basic experimental setup is 
achieved by rolling a two-dimensional topological system 
into a cylinder and threading it with a magnetic flux. 
For our one-dimensional photonic simulator, if we regard 
the left and right edges of the photonic lattice as the two 
ends of the cylinder, the periodic parameter 9 as the ex¬ 
ternal magnetic flux and the in-gap energy as the fermi 
level, our system can be naturally used to simulate the 
experimental setup in Laughlin’s pumping argument and 
to test scattering theory of topological insulators. When 


the frequency of the incident photon towards one edge 
is tuned into the in-gap energy and the external peri¬ 
odic parameter 0 is tuned over one period, the pumping 
particle number per cycle can be expressed as 

where r{0) is the reflection coefficient of the incident pho¬ 
ton from one edge. In this way, based on scattering the¬ 
ory of topological insulators [^, the topological invari¬ 
ant can be characterized b y th e winding number of the 
reflection coefficient phase [ 2 ^. 

To further demonstrate this point, we use Green func¬ 
tion to analytically derive the reflection coefficient from 
the left edge of the above one-dimensional lattice (see 
supplemental materials), giving 


r{e) 


mi + im2 
mi — im2 ’ 


( 7 ) 


where mi = 45J + {Ep + JeCOs{6)){Ac + Jecos(0)) — 
2JJeSm{0) - ^[Ej - 4J2 - J| cos2(6»)][j;2 - 452 - J2 

-|-4Je(5sin(6l)], m 2 = {Ep + Jecos(0))-^J2 - {Ep + Ac)^, 
Ep is the in-gap energy. By substituting the above equa¬ 
tion into Eq. ®, we get 


Q = 


( 8 ) 


One finds that the winding number of the phase of the 
reflection coefficients is exactly equal to the topological 
invariant of this system. In the following section, we 
will show that the information regarding the photonic re¬ 
flection coefficient can be probed spectroscopically using 
cavity input-output process. The photonic Chern insula¬ 
tor is then detected by counting the winding number of 
reflection coefficient phase. 

Probing edge states and reflection coefficients. In con¬ 
trast to fermi system, one can directly probe the edge 
state and its scattering feature in our photonic simu¬ 
lator. The reason is that bosonic photons can occupy 
one particular eigenstate at the same time. This could 
be done by externally driving the resonators with the 
driven frequency tuned as the eigenenergy of the lat¬ 
tice, then the corresponding eigenmode would be oc¬ 
cupied with some weights. In the rotating frame with 
respect to the driving frequency, the driven Hamilto¬ 
nian is Hd = J2ni^naaf + ^nbbf) + h.C, where Vtna.nb 
are the driven amplitudes in the nth unit cell. In 
the presence of dissipation, the expectation value of 
the cavity field aj in steady state can be derived from 
the solution of the Lindblad master equation (dj) = 
—i{[aj,H + Hd]) + Kwhere the Lindblad 
term L[an]aj = anajaf — {afan, k is the cavity de¬ 

cay rate. In the new bases (z = ((oi), ( 61 ),..., (a„), (Ii„))^’’ 
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and ft = {ftia,flib,---,ftna,f'nb)'^^ with Tr representing 
the transposition of matrix, based on the condition of 
the steady state solution {dj) = 0, we can write the ex¬ 
pectation value of the cavity fields in the steady state 
as 


d=-{A, + T-i'^)-% (9) 

where the elements of matrix T are defined by Tna,nb = 
Tnb,na — Jl deSin(0)/2, Tna,{n — l)b ~ ^(ri —l)b,na = <^2 + 
de sin(d)/2, TjiQ^(^b) ,na{b) ~ i</eCOs(^). 


(a) (b) 



- K = 0.1 Je 

-K = O.TJe 

- K.= l.SJe , 



FIG. 3: (a-b) The average photon number in the steady state 
when the left edge state (circle) and one bulk state (square) 
are driven to be occupied by tuning Ac = 0.5Je and 1.2Je, 
respectively. For edge state case, the average photon numbers 
are also plotted when the driven pulse is applied on the middle 
(diamond) and rightmost (star) resonators. The reflection 
coefficients from the left edge for topological (c) nontrivial 
and (d) trivial insulators. The lattice size is (a,c,d) L = 10 
and (b) L = 4, the driving amplitude flia = O.lJe and the 
cavity decay rates k = O.lJe (solid), 0.7Je (dash) and 1.5Je ( 
dash-dot). 

To probe the edge states, we need to occupy this edge 
states firstly. As shown in Fig. 2(a), there is one pair of 
edge states at the in-gap energy for photonic topological 
insulator. In particular, we choose to excite the left edge 
state through external driving the leftmost resonator (see 
Fig. 1(a)), with the driven microwave pulses chosen as 
ft = (riia, 0,..., 0,0)^’' and driven frequency ojd tuned 
to the in-gap energy. The reason is that the left edge 
state has maximal probability occupying the leftmost res¬ 
onator. In Fig. 3(a), we have plotted the average photon 
number in the steady state for this case and find that 
the photons are most localized in the left edge resonator. 
In contrast, if the middle and rightmost resonators are 
driven with same laser, the occupied probability of the 


left edge mode is very small, then there will almost be 
no resonant eigenmode and all the photons will finally 
decay into vacuum in the steady state. In contrast, when 
the driven frequency is tuned as the bulk energy, the 
photons are extensively populated in the lattice, which 
satisfies the feature of Bloch bulk state. The results stay 
even if we choose two unit cells (see Fig. 3(b)). There¬ 
fore, the photonic edge state can be directly observed by 
measuring the corresponding average photon number in 
the steady state. 

The detection of photonic reflection coefficient is nat¬ 
urally related to cavity input-output process [s^. Us¬ 
ing input-output formalism, the reflected output photons 
Qout igf^ edge resonator is related to the input 

photon through a°“* = a™ -I- , where the input field 

g)” is related to the external driving by = iflia 

[3^ . Using Eq. ([9]), the photonic reflection coefficient 
from the left edge is obtained as 

l„out\ 

= 1 + *«,[(A, + T- z-)-i]n. (10) 

\®i / ^ 


In Fig. 3(c) and (d), we plot the numerical results of 
reflection coefficients for photonic topological nontrivial 
(Chern number (7=1) and trivial insulator (Chern num¬ 
ber (7 = 0). The results show that the winding number 
of the reflection coefficient phase of is 1 and 0 respec¬ 
tively, which yield the photonic topological invariants. 
This method also applies for the right edge case and the 
conclusion is same. 

Moreover, we also calculate the corresponding winding 
numbers even when the lattice size is L = 4 (two unit 
cells) and find that the corresponding trajectories are 
same and all the results stay, which means that topo¬ 
logical states could be implemented with only a few- 
resonator lattice. Such remarkable feature is quite attrac¬ 
tive to circuit-QED experimenters. The current circuit 
QED experiment has already realized the capacitive and 
qubit-assisted coupling between two resonators 0 (one 
unit cell). Compared with previous works, our scheme is 
very promising and it is expected to be experimentally 
demonstrated. In all cases, we also take into account the 
influence of the cavity decay. The results show that, if 
the cavity decay rate is not larger than the energy gap 
2 Je, the in-gap energy will remain in the energy gap and 
the winding number will remain the same, then our mea¬ 
surement is very robust to fluctuations of the frequency 
of the input photon. 

Experimental discussion. Before concluding, a detailed 
estimate of the experimental parameters involved is in 
order. For circuit QED experiment [^, with a typical 
choice of ojd = 5A, = O.IA, the qubit-assisted hop¬ 

ping rate Je can approach the order of 10 MHz. For 
the current coupled transmission line resonator experi¬ 
ment l37| . the hopping rate Ji ,2 can be tuned within the 


range 1 — 100 MHz. One can easily check that the exper¬ 
imental parameters required in our work are within the 
experimentally accessible regimes. For the experimen¬ 
tal detection of the reflection coefficient phase, we as- 
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sume an initial input driven lasers on the leftmost (right¬ 
most) resonator prepared in a coherent state |a), then 
the reflected output photon pulse will be in the coherent 
state where (j) is the reflection coefficient 

phase. Based on homodyne detection by interfering the 
output photons with a local oscillator, the phase (j) in 
the output coherent state can be extracted. 

In this way, the winding number of the photonic reflec¬ 
tion coefficient phase is measured, yielding the photonic 
topological invariant. 

Summary. In summary, we have introduced a concep¬ 
tually simple method to realize a photonic Chern insu¬ 
lator in a one-dimensional circuit-QED lattice. Based 
on Laughlins pumping argument and input-output for¬ 
malism, we have further demonstrated that the photonic 
edge states and topological invariant can be unambigu¬ 
ously measured even in a dissipative few-resonator net¬ 


work, which may take a significant step towards observ¬ 
ing topological invariant with circuit-QED. Our method 
also provides a new route and simple means to study and 
probe high dimensional photonic topological states. By 
introducing effective photon interaction, the exotic frac¬ 
tional photonic Chern insulator could be further studied 
in such framework. 
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Appendix A: Tunable flux-qubit-resonator coupling 


<- 



FIG. 4: Setup for tunable flux-qubit-resonator coupling. 

In each unit cell of our lattice system, the setup consists of two resonators coupled by one flux qubit. The flux 
qubit is made of three Josephson junctions. Using the method designed in Hi, the coupling between resonators 
and flux qubit can be mediated by Josephson junctions. In the following, we will focus on the dominant energies of 
Josephson junctions V{4>n) = cos(</>„). Aussume the seven Josephson junction energies are Eji = Ej 2 = Ej, 

Ej 3 = aEj, Ej 4 = Ej 5 = f3Ej and Ejq = Ejr = (3Ej, each loop is applied by an external magnetic field with the 
flux fj = $j/<I>o- Then the total potential energy of the whole system can be written as 

7 

[/ = - ^ Ejn COs{(j)n) 

71 — 1 

Based on the flux quantization condition 

4>i — 4>2 + 4>3 = fi 
^'4’a = /2 
<^5 ~ <^4 = fs 

<1)2 - 4>i + 4>6 - ^i’b = fi 
(i>7 — = /s 


(Al) 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 


we can further rewrite the above potential as 


U = -£’j[cos((/)i) + cos((/)2) + q:cos(/i + (^2 - <;^i) 

+ 213 COs(/ 3 ) COs(/a + 02 - 01 + A0a) 

+ 27 COs(/ 5 ) COs(/b - 02 + 01 - Al/^b)] 


(A7) 

(AS) 

(A9) 
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where fa = fi + f 2 + fs/^, fb = fi + / 5 / 2 , the resonator fields are = Atpo{a + a+) and Ailjb = Atpo{b + b'^). 

By considering the resonator part and omitting the second order qubit-resonator coupling [I|, , based on the two 

degree of freedom (c^i, one can get the following qubit-resonator Hamiltonian 

H = HloqcTz + huJc{a~^a -I- b'^b) + E [SaC-afi {a + a~'') + gbCbfj.{b + b~'')](Tfj, (AlO) 

fi—x,y,z 

where ga = 2l3Ej cos{f 3 )AipQ and gb = 2j3Ej cos{fb)A'ipQ^ Ca^ and Cb^ are determined by the parameters 
(a,/?,/i,/ 2 ,/ 4 ). Similar to the method in Q, the next step is to numerically find the parameters to make sure 


Cay,az = Cby,bz = 0 and Cax,bx = 1, then one can get 

E[ = hwoaz + fujJc{a~^a + b'^b) + [ga{a + a“'') -I- gb{h -I- b'^)\ax (AH) 

Note that the above Hamiltonian becomes the Jaynes-Cummings model after making rotating wave approximations. 
To further get the coupling modulations, the systems needs to work at 

h = 0/2 (A12) 

h = (7r-0)/2 (A13) 

fb + fa = '^/2 (A14) 

then we get the required periodic modulation on qubit-resonator couplings as 

gi=ga = go sin(6l/2), 32 = gb = go cos(6l/2) (A15) 


where go = 2f3EjAijjo. Actually, it is not necessary to have tunable qubit-resonator coupling strengths exactly with 
cosine and sine function forms. Instead, for the same purpose, one only needs to separately tune gi and g 2 - Then, 
by introducing a mixing angle 9 = 2 arctan(gi/g 2 ), 52 and gi can be written in 6 space with cosine and sine forms as 
shown in the above equation. 


Appendix B: Reflection coefficient 


In section, we give details using Green function to derive the reflection coefficient Eq. (10) in the main text. For 
this purpose, we model the total system in the scattering process into three parts: left lead, device (our system) and 
right lead. In particular, the reflection coefficient from the left lead is calculated below. The total lattice Hamiltonian 
of the system is H = Eli, + + Hp -t- TIrd + 77r, where 

L/2 j j L/2 

Hd = - Y sin6>)a)j6„ -b (J 2 + y sm6»)6),a„+i -b H.c.] + ^[Je cos6l(a)ja„ - 6)j6„) -b Ac(a)ja„ -b b^bn)], 

n—1 n—1 

J J 

Hl = - w X! (EiCi + H.c.),i 7 R = -^ E (EiCi + H.C.). 

2^-1 i^L + 1 

(Bl) 


We assume the lattice sites of device L is even. The tunnelings between leads and device are given by 

^LD = -^(ajco + H.C.), iJRD = -^(cj^+i^L + H.C.). (B2) 

In the basis {• • •, cLj^, cj, a\,b\,- ■ -, 0 ^/ 2 ’ ^1/2’ formulate the Hamiltonian of the whole system 

as 


Hl 

tl 

0 ■ 

rl 

Hb 

Tn 

0 

4 

Hr_ 


(B3) 


H = 




where 


Hl = 


tl = 


ffo = 


J^ — 


Sint 


Ac 

J 2 


— Jg cos 9 J 2 + ^ sin 0 


A sint 


2 ^ 

Ac + Je COS 9 
Jt — sin 9 


0 

Ji — ^ sin 9 
Ac — Jp cos 9 


". 





r 0 


-J/2 

0 

. . ■ 

... 0 

-J/2 

0 



-J/2 

0 

-J/2 


••• -J/2 

0 - 

-J/2 

5 

-Wr — 

0 


-J/2 

0 


L • • • 0 

-JI2 

0 J 

00 X 00 











r 0 

0 

0 




0 0 

■ ■ 0 









0 0 

• • 0 


VR = 

0 

0 

0 


1 


-J/2 0 

• • 0 J 

00 X L 


L - J/2 

0 

0 


Lxoo 


Ac -1- Je COS 0 Ji - 

^ sin 0 

0 



0 




J LxL 


(B4) 


Then the Green function for the device is given by 




Gd = [El-Hj,-El-E^^]-\ 


(B5) 


where the self-energies of the leads are E£ = and and the lead Green functions are 

gl = [{E + ig)l-H^]-\ 
g^^ = [{E + zg)l-H^]-\ 


(B6) 


After some straightforward calculations, we find that the non-zero elements in the self-energies are [E l]ii = "tISlIoo.oo 
and = ^[g^ii-, otherwise is zero. Due to the symmetrical configuration of the whole system, we note that 

[E£]ii = [E£^]Li. Furthermore, the dispersion relation of the semi-infinite lead is i? = — Jcosfc, the group velocity in 
the lead is = §f = Jsinfc and the self-energy of the lead is [E£]ii = [Ej£]ii = — 0. Keeping vl > 0 

and vii < 0 for the photon injecting from leads to device, we have (assume J > 0 from now on) 


- E^, El = [E£]ii = 1{E- WJ^-E^), 

Er = [E^]ll = ^{E + ly'j^-E^). 

Therefore, the device Green function is 


(B7) 


Gd 


E — Ac — Je COS 9 — 
—Ji -I- ^ sin0 
0 
0 
0 


— Ji -I- ^ sind 0 

E — Ac + Je COS 9 — J 2 — ^ sin 0 
- J 2 - A sin 0 
0 

0 0 


0 

0 


E — Ac — Je COS 0 
—Ji -I- A sin 0 E 


0 

0 

0 

— Ji -I- A sin 0 
Ac -I- Je cos 0 — Er 


-1 


LxL 

(B8) 


Via the continued fraction method and taking into account the periodic pattern of the matrix elements in Gd, the 
closed form of [Gd]ii can be obtained, 


[GdIi/ 4” '^L — E — Ac — Je cos 0 — 


(Ji — ^ sin0)^ 


E — Ac + Je cos 0 


iJ 2 + ^ sin 6)^ 


(B9) 


Solving this algebra equation we have 


[Gd]ii = — 


2{Ep + Pi) 

mi — 17712 ’ 


(BIO) 
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where 


mi = Jf - J| + {Ep + Pi)(A, + Pi) - (Ji + J2)P2 - 

m 2 = {Ep + Pi)^J2-(Pp + A,)^ 

Pi = Je COS 9, P 2 = Je sin 9, 


(P 2 - Jl + J2)2 - p2], 

(Bll) 


and Ep = E — Ac is the in-gap energy of our photonic system. Based on the Fisher-Lee relation Q, Snm{E) = 


—^nm + iy/EnJ^lGujnm, where the scattering matrix is S' = 
thus 


J'L iR 
th rn 


. The reflection coefficient from the left lead is 


rL{9) = -l + tVj^-E^Gu]ii. 

Therefore, through substituting Eq. (IBIOI) into Eq. (IB12L we find that 


rL{9) = 


mi + im2 

mi — im2 


(B12) 


(B13) 
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